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ABSTRACT 



Aims. Based on the archival data from the Chandra observations of nearby galaxies, we study different sub populations of low- 
mass X-ray binaries (LMXBs) - dynamically formed systems in globular clusters (GCs) and in the nucleus of M31 and (presumably 
primordial) X-ray binaries in the fields of galaxies. Our aim is to produce accurate luminosity distributions of X-ray binaries in 
different environments, suitable for quantitative comparison with each other and with the output of population synthesis calculations. 
Methods. Our sample includes seven nearby galaxies (M31, Mafi"ei 1, Centaurus A, M81, NGC 3379, NGC 4697, and NGC 4278) and 
the Milky Way, which together provide relatively uniform coverage down to the luminosity limit of 10^^' erg/s. In total we have detected 
185 LMXBs associated with GCs, 35 X-ray sources in the nucleus of M31, and 998 field sources of which ~ 365 are expected to be 
background AGN. We combine these data, taking special care to accurately account for X-ray and optical incompleteness corrections 
and the removal of the contamination from the cosmic X-ray background sources, to produce luminosity distributions of X-ray binaries 
in different environments to far greater accuracy than has been obtained previously. 

Results. We found that luminosity distributions of GC and field LMXBs differ throughout the entire luminosity range, the fraction of 
faint (log(Lx) < 37) sources among the former being ~ 4 times less than in the field population. The X-ray luminosity function (XLF) 
of sources in the nucleus of M31 is similar to that of GC sources at the faint end but differs at the bright end, with the M31 nucleus 
hosting significantly fewer bright sources. We discuss the possible origin and potential implications of these results. 

Key words. X-rays: binaries ~ (Galaxy:) globular clusters: general 



1. Introduction 

It has long been known that there are many more low-mass X- 
ray binaries (LMXB^ per unit stell ar mass in Galactic globular 
clusters (GCs) than in the field (Clark !T975l) . This fact is conven- 
tionally explained as a result of dynamical formation of LMXBs 
in the high stellar density environment of GCs where the prob- 
ability of two-bod y interactions, which scales as pi, is high 
dFabian et al.ll975l) . In the Chandra era this picture received fur- 
ther support from the high specific frequenc y of LMXBs in GCs 
observed in nearby e xternal galaxi es (e.g., Aneelini et al. 2001; 
ISarazin et al.] 120031; iMinniti et al.] 12004; Jordan et al.. .2007b) . 
Also, a significant "surplus" of LMXBs was detected in the nu- 
cleus of M31, with the spatial distribution of compact X-ray 
sources following the "p^" law (Voss & Gilfanov 2007a). The 
stellar density is low outside of GCs and the nuclear region of 
galaxies, with a correspondingly lower probability of stellar in- 
teraction, therefore primordial formation is thought to be the 
main formation process for LMXBs in the main bodies of galax- 
ies. Their volum e densities follow the distribution of stellar mass 
(IGilfanovl[200l) . 



' Throughout this paper we refer to objects that have been actively 
accreting in recent times (i.e. \og{Lx) i 35) as X-ray binaries. 



Although the above picture is attractive in its simplic- 
ity, there is a plausible alternative scenario: The entire pop- 
ulation of LMXBs in galaxies, including those in the field 
may have been produced dynamically in GCs and later ex- 
pell ed into the field. Although the debate is still going on 
(W hite et al.ll2002t iKundu et alJl2002L 120071: 1 IrwinI 120051: Lluettl 
2005 : [Humphrey & Buotd 20081) . several strong arguments have 
been presented which suggest a (significant) fraction of field 
LMXBs formed in situ via primordial binary formation. These 
include t he difference in spat ial distributions of field LMXBs and 
GCs (e.g. lKundu et"ani2007h and the lack of correlation between 
the specific frequency of field LMXBs and that of GCs (e.g. 
Juett 2005). At the same time, Humphrey & Buote (2008) came 
to the opposite conclusion. The recently found evidence that 
the X-ray luminosity functi on (X LF) of GC - LMXBs may differ 
from that of field LMXBs & Gilfanovl [20075 IVoss et al.1 
2009:'Woodle y et al.ll2008t iKim et al.i 20091 adds to this debate. 
Although some caveats are in order, differences in the luminos- 
ity distributions of the GC and field binaries suggest that the two 
sub populations of LMXBs may hav e different formation and/or 
evolution histories (IVoss et al.ll20"09h . 

Differences in the luminosity distributions of LMXBs in GCs 
and in the field may be most obvious in the low-luminosity 
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Table 1. The sample of external galaxies observed by Chandra 



Galaxy 


Type 


Distance 


N„ 


Study Field 




M, 


Exp 


Sensitivity 


dx/dy 


Conversion Factor 






(Mpc) 






(Mo/Lk,o) 


(lO'^Mo) 


(ks) 


(erg s"') 


(pixel) 


(erg cm"^ count"') 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


M31 


Sab 


0.78+0.03 


6.7 


r = 11' 


0.56 


2.7 


-300 


4x10^'" 


+0.02/-0.37 


3.4 X 10"^ 










15', 9', -65° 




0.60 


-150 


7x10^* 


+0.11/+0.51 




Cen A 


SO 


3.4±0.4 


8.6 


r = 10' 


0.76 


6.4 


-800 


6x10" 




3.5 X 10-' 


M81 


Sab 


3.63±0.34 


4.2 


10', 5', -20° 


0.70 


6.0 


-240 


7x10" 


+0.10/+0.49 


3.2 X 10-' 


Maffei 1 


SO 


3.0±0.3 


85.1 


D25 & HST 


0.73 


1.1 


-55 


3x10^6 


+0.01/+0.29 


7.7 X 10-' 


N3379 


El 


11.1 


2.8 


D25 & HST 


0.83 


5.4 


-330 


4x10^'^ 


-h 


3.1 X 10-' 


N4697 


E6 


11.8 


2.1 


D25 


0.77 


5.8 


-200 


5x10^*^ 


-/- 


3.0 X 10-' 


N4278 


El-2 


16.1 


1.8 


D25 & HST 


0.79 


4.2 


-480 


6x10^*^ 


-1- 


3.0 X 10-' 



(1) - Galaxy name. For M31, the first line is for the bulge region, the second line is for the region in the disk. (2) - Galaxy Type. (3) - Distanc e 
and its uncertainty (when available). References and methods are: M31 - luminosity function of red clump stars (IStanek & Gamavich"1998|); 
Centaurus A ~ Cepheids (Ferrarese et al. 2007); M81 - Cepheids ([pr eedman et al. 1994); Maffei 1 - galaxy fundamental plane (Fingerhut et al. 
2003{) ; NGC 3379 - luminosity function of GCs ( Kundu & Wh itmoreli2001l) ; NGC 4697, NGC 4278 - surface brightness fluctuation (iTonrv et al. 
2001). (4) - Galactic column density (Dickev & Lockman 1990). (5) - The region used to study XLFs. When three numbers are given, they refer 
to major, minor axi s and position a ngle. (6) - A^-band mass-to-light ratios derived from Bell & de Jong (2001), with B -V colors from the RC3 
catalog ( de Vaucouleurs et al, ||1991|) except for Mafl'ei 1, which is fromButa & McCall ( 1983). (7) - Stellar mass in the study field, as calculated 
from the A^-band magnitudes derived from 2MASS Large Galaxy Atla s (jjarrett et al. 2003). For M31 we used the IRAC/Spitzer data, and the 3.6 
ytim flux was converted to /f-band foUowing lBogdan & Gilfanovl ( l2010h . (8) - The total exposure time of Chandra observations. (9) - Point source 
detection sensitivity estimated from the incompleteness functions in Fig.jS] (10) - Attitude correction. (1 1) - Conversion factor of Chandra count 
rate to unabsorbed X-ray flux in the 0.5-8 keV band. 



(log(Lx) - 37) domain. Thus the reliable detection and quan- 
titative study of any possible differences in the XLF demands 
special care in the treatment of incompleteness effects and the re- 
moval of the cosmic X-ray background sources (CXB). Another 
difficulty, of a more fundamental nature, is the statistical noise 
caused by the small numbers of sources. Although the majority 
of previous investigations s eem to converge i n their conclusions, 
with a few exceptions (e.g.. lVoss et al.ll2009h most of these stud- 
ies have marginal statistical significance. However it is difficult 
to achieve higher quality statistics by studying individual galax- 
ies that host a limited number of sources. Massive ellipticals with 
their large GC populations could avoid this difficulty. But the 
long distances to the best candidates require deep X-ray obser- 
vations, in the Msec range, to reach the required depth. Such 
data sets are not available in the Chandra archive. However, the 
potential impact of accurately determining the luminosity distri- 
butions of X-ray binaries located in different environments on 
our understanding of the formation and evolution of LMXBs is 
high. This motivated us to attempt to produce the most accu- 
rate LF of GC-LMXBs to date by combining Chandra data for 
multiple galaxies. To this end we undertook a systematic survey 
of nearby galaxies with sufficient numbers of LMXBs and GCs. 
The results of this study are reported below. The paper is struc- 
tured as follows. We describe our selection criteria and resulting 
sample in section|2]and the data preparation and analysis in sec- 
tion [3] In section |4] we describe the procedure used to combine 
XLFs. Results are presented and analyzed in section |5] and dis- 
cussed in section|6l Section|7]lists the conclusions. 

2. The sample 

In constructing the sample our goal is to provide uniform cov- 
erage over as wide a range in luminosity as possible. We aim to 
study sources as faint as lO''^ erg/s. On the other hand, our goal 
is to have good enough statistics at the bright end where the spe- 
cific frequencies of sources (per GC or per unit stellar mass) is 
low. Therefore our strategy is to include all galaxies with the best 



sensitivities achieved by Chandra so far and complement this 
with several sufficiently massive galaxies with somewhat low 
sensitivity in order to properly sample the high-luminosity do- 
main. We based our selection on the list of normal galaxies avail- 
able in the public Chandra archive. We did not exclude late-type 
galaxies, but in constructing the XLF of the field sources we con- 
sidered only their bulges (to exclude possible contamination by 
HMXBs). The main selection criterion used was a detection sen- 
sitivity better than log(Lx) ~ 36.5-37. This translates into a joint 
constraint on the distance to the galaxy and the exposure time 
of the Chandra observation. We also decided to exclude galax- 
ies with stellar mass less than lO'^M© because of their smaller 
LMXB populations and the consequently higher contamination 
by resolved CXB sources. Finally, we required the availability 
of extensive GC data in order to reliably separate GC and field 
sources. Our final sample includes seven nearby galaxies (Table 
[T]i- In addition we also include the GC sources in the Milky Way. 

M3 1 is the only nearby galaxy with Chan dra sensitivity bet- 
ter than 10^^ erg s-'. IVoss & Gilfanovl (12007 al) analysed 160ksec 
of Chandra data available at the time and found 12 LMXBs in 
confirmed GCs in the bulge. Since this study, an additional - 140 
ksec of data has been collected by Chandra (which brings the 
total exposure time of the bulge to ove r 300 ksec), and mor e ac- 
curate GC data have been published (iPeacock et al.ll20T0 ). We 
also analysed an additional 160 ks observation of a region in 
the disk. Centaurus A was the target of a recent Chandra VLP 
program. With a total Chandra exposure time of -800 ks, a de- 
tection sensitivity of 6x1 0^^ erg s-' has been reached in this 
galaxy. ' Voss et alJ (12009^ find 47 GC-LMXBs in this galaxy, 
so we use their source lists in our analysis. A similar detection 
sensitivity was reached in M8 1 with an exposure time of -240 
ks. The four other external galaxies in our sample have detec- 
tion sensitivities of a few x 10^^ erg s-' and are included to in- 
crease the statistics of luminous sources. One of them, Maffei 
1, is relatively small and marginally passed our mass thresh- 
old. However, it appears to be particularly rich in X-ray sources. 
The X-ray populations in NGC 3379, NGC 4697 and NGC 4278 
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Table 2. The list of Chandra observations analyzed in this paper. 



Galaxy 


Obs-ID 


Instrument 


Exp. (ks) 


Galaxy 


Obs-ID 


Instrument 


Exp. (ks) 


Galaxy 


Obs-ID 


Instrument 


Exp. (ks) 


M31(l) 


0303 


ACIS-I 


12.01 


M31(l) 


7064 


ACIS-I 


29.07 


M81 


5937 


ACIS-S 


12.16 


M31(l) 


0305 


ACIS-I 


4.18 


M31(l) 


7068 


ACIS-I 


9.62 


M81 


5938 


ACIS-S 


11.96 


M31(l) 


0306 


ACIS-I 


4.18 


M31(l) 


7136 


ACIS-I 


4.96 


M81 


5939 


ACIS-S 


11.96 


M31(l) 


0307 


ACIS-I 


4.17 


M31(l) 


7137 


ACIS-I 


4.91 


M81 


5940 


ACIS-S 


12.13 


M31(l) 


0308 


ACIS-I 


4.06 


M31(l) 


7138 


ACIS-I 


5.11 


M81 


5941 


ACIS-S 


11.96 


M31(l) 


0309 


ACIS-S 


5.16 


M31(l) 


7139 


ACIS-I 


4.96 


M81 


5942 


ACIS-S 


12.11 


M31(l) 


0310 


ACIS-S 


5.14 


M31(l) 


7140 


ACIS-I 


5.12 


M81 


5943 


ACIS-S 


12.17 


M31(l) 


0311 


ACIS-I 


4.96 


M31(l) 


8183 


ACIS-I 


4.95 


M81 


5944 


ACIS-S 


11.96 


M31(l) 


0312 


ACIS-I 


4.73 


M31(l) 


8184 


ACIS-I 


5.18 


M81 


5945 


ACIS-S 


11.72 


M31(l) 


*1575 


ACIS-S 


38.15 


M31(l) 


8185 


ACIS-I 


4.95 


M81 


5946 


ACIS-S 


12.17 


M31(l) 


1577 


ACIS-I 


4.98 


M31(l) 


8191 


ACIS-I 


4.95 


M81 


5947 


ACIS-S 


10.84 


M31(l) 


1581 


ACIS-I 


4.46 


M31(l) 


8192 


ACIS-I 


5.09 


M81 


5948 


ACIS-S 


12.18 


M31(l) 


1582 


ACIS-I 


4.36 


M31(l) 


8193 


ACIS-I 


5.16 


M81 


5949 


ACIS-S 


12.18 


M31(l) 


1583 


ACIS-I 


5.00 


M31(l) 


8194 


ACIS-I 


5.04 


M81 


9122 


ACIS-S 


10.04 


M31(l) 


1585 


ACIS-I 


4.95 


M31(l) 


8195 


ACIS-I 


4.95 


Maffei 1 


5619 


ACIS-S 


55.75 


M31(l) 


1854 


ACIS-S 


4.75 


M31(2) 


0313 


ACIS-S 


6.05 


N3379 


1587 


ACIS-S 


31.92 


M31(l) 


2895 


ACIS-I 


4.94 


M31(2) 


0314 


ACIS-S 


5.15 


N3379 


*7073 


ACIS-S 


85.18 


M31(l) 


2896 


ACIS-I 


4.97 


M31(2) 


1576 


ACIS-I 


4.95 


N3379 


7074 


ACIS-S 


69.95 


M31(l) 


2897 


ACIS-I 


4.97 


M31(2) 


1580 


ACIS-S 


5.13 


N3379 


7075 


ACIS-S 


84.18 


M31(l) 


2898 


ACIS-I 


4.96 


M31(2) 


1584 


ACIS-I 


4.97 


N3379 


7076 


ACIS-S 


70.14 


M31(l) 


4360 


ACIS-I 


4.97 


M31(2) 


2049 


ACIS-S 


14.76 


N4697 


784 


ACIS-S 


39.76 


M31(l) 


4678 


ACIS-I 


4.87 


M31(2) 


2050 


ACIS-S 


13.21 


N4697 


*4727 


ACIS-S 


40.45 


A>n 1 / 1 \ 


A Aid 




"t.l 1 




ZUjl 




1 T on 


iS4Dy / 


A '70Q 
4/ZS 




1 ^ 

ib. lo 


M31(l) 


4680 


ACIS-I 


5.24 


M31(2) 


2894 


ACIS-I 


4.72 


N4697 


4729 


ACIS-S 


38.61 


M31(l) 


4681 


ACIS-I 


5.13 


M31(2) 


2899 


ACIS-I 


4.97 


N4697 


4730 


ACIS-S 


40.58 


M31(l) 


4682 


ACIS-I 


4.93 


M31(2) 


2901 


ACIS-I 


4.68 


N4278 


4741 


ACIS-S 


37.94 


M31(l) 


4719 


ACIS-I 


5.10 


M31(2) 


2902 


ACIS-I 


4.76 


N4278 


'1011 


ACIS-S 


111.72 


M31(l) 


4720 


ACIS-I 


5.14 


M31(2) 


•4536 


ACIS-S 


54.94 


N4278 


lOlS 


ACIS-S 


52.09 


M31(l) 


4721 


ACIS-I 


5.16 


M81 


•0735 


ACIS-S 


50.56 


N4278 


7079 


ACIS-S 


106.42 


M31(l) 


4722 


ACIS-I 


4.87 


M81 


5935 


ACIS-S 


11.12 


N4278 


7080 


ACIS-S 


56.54 


M31(l) 


4723 


ACIS-I 


5.05 


M81 


5936 


ACIS-S 


11.55 


N4278 


7081 


ACIS-S 


112.14 



M31(l) is the bulge region and M31(2) is the disk region. The observations marked by "*" were used as the reference when combining the data. 



have been studied previously (e.g. jKundu et al.ll2007tlKim et al.l 
[2009; Brassington et al. 2009). For these galaxies we redid the 
data analysis and found it to be in overall agreement with the 

above authors. 

The Milky Way hosts 150 GCs (lHarrislll996h . of which 12 
are known to host bright LMXBs. As all the Milky Way GCs 
have been surveyed in the X-rays multiple times by various in- 
struments, we assume that our sample of GC-LMXBs is com- 
plete. We used the data from the All-Sky Monitor aboard RXTE 
to measure the luminosities of these sources. The advantages and 
shortcomings of such an approach are discussed in section U4l 

3. Data analysis 

3. 1 . Source detection 

The details of the Chandra observations of external galaxies 
used here are listed in Table |2] These observations were re- 
duced following the standard CIAO threads (CIAO version 3.4; 
CALDB version 3.4.1). Time intervals when background flares 
happen were not excluded, since the benefit of the increased ex- 
posure time outweighs the increased background. The energy 
range was limited to 0.5-8.0 keV. We have only used sources lo- 
cated in the regions where Chandra data overlaps with the GC 
data to construct the XLFs. The sizes of these "study fields" are 
listed in Table [1] and the regions are overlaid on the X-ray im- 
ages shown in Fig.[T] 

In order to detect sources we used "wavdet ect" with the same 
parameters as IVoss & Gilfaiiovl (l2006i l2007ah . Thresholds were 



set to be 10"^ for all galaxies, yielding an average of 1 false 
source per 10^ 0.492x0.492 arcsec^ ACIS pixels. In order to 
prevent distortion of point sources due to attitude reconstruction 
errors in the course of combining images, we used the bright- 
est sources detected within a 4' radius of the telescope axis in 
each observation and corrected offsets between observations us- 
ing the method described in Voss & Gilfanov (200 7a). For the 
disk region of M31, the observations are distributed in too wide 
an area to make this correction possible. For M81 and NGC 
4697, the offsets between observations were insignificant thus 
this step was ski pped. For NGC 3379, we used offsets from 
iBrassington et al.l (1200 8). The observations were then shifted to 
match the coordinate system of the reference (marked with an 
asterisk in Table |2]l and thereafter combined together and re- 
analyzed. 

The next step is to apply an absolute astrometry correction to 
the combined image. We used the 2MASS All-Sky Point Source 
catalog ( Skrutskie et al. 2006) and correlated it with the brightest 
X-ray point sources. X-ray images were shifted to give the short- 
est rms-distances between the X-ray sources and their 2MASS 
counterparts. These corrections are listed in Table [1] where dx 
refers to correction in the west and dy is the correction to the 
north. In the case of three galaxies: NGC 3379, NGC 4697, and 
NGC 4278, this step was skipped because we did not find enough 
matches. 

The statistics of the detected sources are summarized in 
Tables [3] and H] To estimate the count rates we followed the 
method described in IVoss & Gilfaiiovl (I2007ah . The luminosities 
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Fig. 1. X-ray images (0.5-8 keV) of external galaxies and the fields of view studied. For M31 and M81 we study the sources within 
the ellipses with radii and position angle shown in Table[T] For Maff'ei 1, NGC 3379, NGC 4697, and NGC 4278 ellipses show D25 
regions and squares show the HST fovs. The field of view of interest are the overlapping regions between the two. Crosses show 
detected GC-LMXBs in each galaxy. The Chandra image of Centaurus A can be found in Fig. 1 in iVoss et al.l (l2009h . 



of point sources were calculated assuming a power-law spectrum 
with r= 1 .7. Count rates in 0.5-8.0 keV band were converted into 
absorption corrected fluxes assuming Galactic absorption. We 
list the conversion factors in Table \T\ 

To estimate the expected numbers of CXB sources among 
detected X-ray sources, w e used the res u lts of the CXB log(A^)- 
log(5') determination by iMoretti et al.l (l2003l) . We used their 
equation (2) for the soft band counts and converted the flux to 
the 0.5-8.0 keV band, assuming a power-law spectrum with a 
photon index of 1 .4. The predicted numbers of CXB sources ac- 
count for the incompleteness are described in section 14.11 The 
computed numbers of CXB sources are listed in Table |3] In the 
closest galaxies, background AGN account for a large fraction 
of detected X-ray sources, especially in M3 1 where nearly half 
of the X-ray sources are CXBs. Maffei 1 is a small galaxy that is 
abundant in LMXBs, and the contamination by CXB sources is 
minimal. In NGC 3379 and NGC 4278 HST the WFPC2 field- 
of-views (FOVs) are located in the very central region, where the 
CXB fraction is less than 10%. In NGC 4697 the CXB fraction 
is about 15% in the D25 region. 

The CXB estimates based on the average source counts are 
subject to uncertainties caused by angular fluctuations of the 
density of background AGN. These are likely to be reduced in 
our analysis as it covers a rather large solid angle composed of 
non contiguous fields. Nevertheless, for each individual galaxy 



we verified that the obser-ved density of compact sources out- 
side its main body is consistent, within the statistical errors, 
with the predicted density of CXB sources. This was possible 
to do direcdy for Maffei 1, NGC 3379, NGC 4697, and NGC 
4278, thanks to their relatively small angular size. For M31 
and Centaurus A whose angular extent excee ds or i s compara- 
ble to th e C handra FOV , we u sed the results o f|\^& Gilfanoi 
(l2007bh and lVoss et al.1 (l2009l) . In CentaurTis A the CXB source 
density was found to exceed the average source count by a factor 
of ~ 1.5, which that was accounted for in our calculations. 

3.2. Identification of GC sources 

We correlated the lists of detected X-ray point sources with the 
GC lists available for M31, M81, Maff^ei 1, NGC 3379, NGC 
4697, and NGC 4278. For Centaurus A we used the results of 
IVossetal.ld2009h . 

For M31, the most recent and complete GC catalog is a sys- 
tematic survey using W FCAM on the United Kingdom Infrared 
Telescope and SDSS bv lPeacock et alJ (l20IOl) . In total there are 
416 confirmed GCs, with 121 located i n the two Chand r a field s 
in our study. GCs in M8 1 are from Perelmuter et al.l (119951) . 
Chandar et al., (1200 ih . and Schroder et al. (2002). Chandar's 
study is based on deep HST observations that cover 25% of our 
Chandra field. We took the 59 GCs from IChandar et al.l (1200 ll) 
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Table 3. Statistics of compact sources I 



Galaxy 


^XRS 




A'gc 




A'gc-x 


"■opt 


R 


A'r 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


M31 


386 


194 


121 


1.00 


26 


1.00 


2.0" 


0.2 


Cen A 


231 


64 


479 


0.67 ± 0.03 


47 


0.55 ±0.19 


2.0" 


1.2 


M81 


220 


79 


77 


0.77 ± 0.09 


8 


0.84 ± 0.48 


3.0" 


0.8 


Maffei 1 


38 


1 


20 


1.00 


4 


1.00 


1.0" 


0.3 


N3379 


59 


4 


61 


0.80 ±0.10 


9 


0.74 ± 0.39 


1.0" 


0.6 


N4697 


117 


17 


441 


0.85 ± 0.04 


39 


0.93 ±0.31 


0.8" 


1.1 


N4278 


120 


6 


266 


0.69 ± 0.04 


40 


0.86 ± 0.29 


0.6" 


1.8 


MW 






150 


1.00 


12 


1.00 






Total 


1171 


365 


1615 




185 






6.7 



Columns are: (1) - Galaxy name. (2) - Total number of resolved X-ray point sources in the study fields. (3) - Predicted number of CXB sources 
in the study fields above the corresponding sensitivity threshold. (4) - Number of optically identified GCs. (5) - Completeness fraction of GC 
lists and its Icr uncertainty, estimated as described in section [33] (6) - Number of LMXBs found in GCs. (7) - Completeness fraction and its 
uncertainty of GC lists with respect to GCs containing LMXBs (see section l33t . (8) - Search radius to match XRS to GC. (9) - Expected number 
of random coincidences of X-ray sources with GCs. 

Table 4. Statistics of compact sources II 



Galaxy 


A'' 

^'gc-x 


^*GC-X 


A'^ 

^'gc-x 




^F-X 


^*F-X 


"f-x 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


M31 


2 


11 


12 


0.10 ±0.03 


110 


64 


28 


9.5 ± 1.8 


Cen A 





16 


30 


0.06 ±0.01 


6 


85 


29 


7.2 ± 1.3 


M81 





4 


4 


0.05 ± 0.03 










Maff'ei 1 





2 


2 


0.10 ±0.07 





6 


12 


14.8 ±4.3 


N3379 








8 


0.13 ±0.05 





9 


24 


6.3 ± 1.3 


N4697 





2 


34 


0.08 ±0.01 





2 


66 


8.9± 1.1 


N4278 





4 


36 


0.14 ±0.02 





3 


52 


17.4 ±2.4 


MW 


1 


9 


2 


0.013 ± 0.009 










Total 


3 


48 


128 


0.08 ± 0.01 


116 


169 


211 


9.8 ±0.7 



Columns are: (1) - Galaxy name. Columns (2)-(4) and (6)~(8) - Number of GC-LMXBs (A'gc-x) and field LMXBs (Wp-x) in different luminosity 
ranges (1,2 and 3 refer to log(Lx) ranges of 35-36, 36-37 and >37) with incompleteness higher than 0.5. The source numbers are not corrected for 
incompleteness, and the CXB contribution is not subtracted. Columns (5) and (9) - The specific number of GC-LMXB (per GC) and field LMXBs 
(per 10'" Mq) in the highest luminosity bin log(Lx) > 37. The numbers are coiTected for incompleteness of X-ray source lists, the contribution 
of CXB is subtracted. Note that the specific numbers of GC-LMXB are not corrected for incompleteness of the GC lists and are given here as a 
characterization of our sample, rather than of the properties of GC systems in different galaxies.). 



and the others from the other two catalogs, which resulted in 77 
confirmed GCs in the Chandra field. For MafFei 1, there are 20 
GCs from HST observation by Buta & Mc Call (2003). F orNGC 
3379, we took the 61 GCs from iKundu & Whitmorg (l200lh . 
which are based on deep HST observations. The GC list for NGC 
4697 is taken from Jordan et al. (201 1). And the GCs in NGC 
4278 a re from lKundu & Whitmord (1200 ih and|Brassingtonetal] 
(l2009l) . 

The search radius R used in cross-correlating the X-ray 
source lists and GC catalogs was chosen for each galaxy indi- 
vidually based on the following considerations. The number of 
random matches is A^i- = jtR^ x A^xrs x Nqc/A, where A^xrs is 
the number of X-ray point sources, A^gc the number of GCs, and 
A the area of our study field. Because of the rather high source 
density, the number of random matches may be non-negligible 
for high values of the search radius R. On the other hand, the 
search radius has to be broad enough to account for position 
errors and, for the closest galaxies, the finite angular sizes of 
GCs. We therefore devised a procedure in which we varied R 
from to 5" (Fig. |2]l. For each value of R we computed the 
number of true matches as the number of total detected matches 
minus the predicted number of false matches calculated from 



the above formula. This number increases with R and saturates 
at some value of R that depends on the typical positional error 
and angular extent of GCs. This value of R may be chosen as 
the optimal match radius. In some cases, however, it results in 
too high a fraction of false matches in the sample. We therefore 
set an additional requirement that the predicted number of false 
matches does not exceed 5% of the total number of matches. 
This procedure is a simplified version of the method used in 
IShtvkovskiv & Gilfanovl(l2005h . The optimal search radius used 
for the program galaxies are listed in Table |3] As expected, there 
is a general trend that nearby galaxies require larger search radii. 

The numbers of X-ray sources associated with GCs are listed 
in Table|3]along with the predicted numbers of false matches. 



3.3. Incompleteness of GC lists 

Although the availability of the high quality GC optical data was 
one of the criteria in selecting our galaxy sample, the GC lists are 
not 100% complete for all of them. The incompleteness of these 
lists can result in incompleteness of the GC-LMXB lists and can 
compromise the shape and (less importantly) the normalization 
of the GC XLF 
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0.5 1 
R (arcsec) 

Fig. 2. Examples of the determination of the optimal search ra- 
dius for cross-correlation of X-ray source lists with GC cata- 
logs for two galaxies - M3 1 (solid lines) and NGC 4278 (dashed 
lines). Upper panel: the number of true matches (Mme, computed 
as the difference between the number of total matches Motai the 
number of expected random matches N^) as a function of the 
search radius R. Lower panel: the ratio of the number of random 
matches to the number of total matches. The dash-dotted line 
is the 5% level. The vertical lines in the upper panel show our 
choice of the search radii for these two galaxies. 




Fig. 3. The observed luminosity functions of GCs for six galax- 
ies in our sample and their best-fit models. The turnover lumi- 
nosity of the model in different bands was fixed at the values 
determined elsewhere (see text), the width cr was a free parame- 
ter of the fit. 



In order to estimate the completeness fraction of the GC lists 
we used the fact that optical luminosity function of GCs (GCLF) 
can be described to good accuracy by log-normal distribution in 
the form 



dN 
dM 



1 



^/27T^ 



■ exp 



(M-ijf 



2cr2 



(1) 



where M is the absolute magnitude of GC, the turnover lu- 
minosity, and cr the dispersion. The turnover luminosity is re- 
markably constant in different galaxies. We used the follow- 
ing values for different bands: fiy - -7.41, //j - -8.46 
(iKundu & Whitmordl200lh and = -7.2 (iJordan et al.ll2007ah . 
The reddening corrected photometry data of GCs for each galaxy 
(see references in section 13.2) was fit by this model using max- 
imum likelihood method. The fitting was performed using only 
GCs above the completeness limit of the optical data for each 
galaxy. The width of the distribution cr and normalizations were 
free parameters of the fit. The data along with the best-fit model 
are shown for six galaxies from our sample in Fig. |3] The com- 
pleteness fraction of the GC lists K^^, was then determined as 
a ratio of the total number of detected GCs (of all magnitudes) 
to the total number predicted by integrating the best-fit model. 
The results are listed in Table |3] Given the completeness limits 
of M3 1 and Maffei 1 data, their GC lists are complete. The list of 
GCs in the Milky Way is also believed to be reasonably complete 
dHarris 2001). In agreement with this, the best-fit values obtained 
for these three galaxies are consistent, within errors, with 1. We 



therefore set the completeness fraction for these three galaxies 
equal to unity. 

The second factor, required to correctly computing GC XLF 



is K 



cc.x 

opt 



the completeness fraction of the GC-LMXB identi- 



fications, caused by the incompleteness of the overall GC lists. 
If the probability of finding an LMXB in a globular was inde- 
pendent of its optical luminosity, the two quantities would co- 
incide: K^p,'^ - Kgpj would hold. However, it has been shown 
that X-ray sources te nd to be associated with brighter GCs (e.g., 
ISivakoffetal.ll2007h . This is illustrated by the Fig. |4] where we 
plot the combined LF of all GCs and all the GCs hosting an X- 
ray source in the three galaxies which GC lists are complete - 
MW, M31 and Maffei 1. To determine Kf,^;^ we assume that LF 
of GCs hosting an X-ray source is the same in all galaxies. Using 
the combined LFs in Fig|4]as the template, we then use the ratio 
of the numbers of GCs hosting X-ray sources above and below 
the threshold magnitude of V = -7 (V = -8 for Cen A) to esti- 
mate the number of missed X-ray sources in GCs in each galaxy. 
The threshold magnitude was chosen so that the GC lists are 
complete above its value. The results of this calculation, along 
with their uncertainties are listed in Table |3] 



3.4. The Milky Way sources 

The luminosities of the Milky Way sources were calculated from 
the light curves measured by the All-Sky Monitor aboard RXTE. 
The light curves were averaged over the period from January 
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0.4 r 



0.3 r 




0.2 r 



0.1 7 



0.0 



Fig. 4. The combined LFs of all GCs in the Milky Way, M3 1 , and 
Maffei 1 (the shaded histogram). The thin solid histogram shows 
the combined LF of GCs hosting X-ray sources. The Poisson 
errors for the latter are indicated by the vertical error bars. The 
statistical errors are ~2-3 times smaller for the combined LF of 
all GCs. 



1996 to June 2009. The count rates were transformed to the 0.5- 
8.0 keV band fluxes with the conversion factor obtained using 
PIMMS: 1 count s"' = 4.3x10""' erg cm"^ s"'. A power-law 
with r=1.7 was assumed. To compute sour ce luminosities we 
used distances to the GCs from'Harrisl (ll996l) . 

The following comments regarding the determination of the 
luminosities of the Milky Way sources are in order. The ASM 
fluxes are averaged over a significantly longer time scale than 
the Chandra data for external galaxies. Although both the ASM 
(~ years) and Chandra {~ 1-10 days of total integration 
time, ~ years time span of observations) integration time scales 
are much longer than the characteristic time scales of the ac- 
cretion disk in these sources the variability of the X-ray light 
curves could in principle result in "clipping" of the XLF i.e. 
smoothing out the extrema of long term variability. The ef- 
fe ct of averaging ASM light curves on the XLF was studied 
bv lPostnov & Kuranovl ( l2005l) . who came to the conclusion that 
flux probability distribution functions for persistent galactic X- 
ray binaries are such that light curve averaging does not modify 
the shape of the power-law luminosity distribution. To verify this 
further we considered variations in the XLF obtained by averag- 
ing ASM light curves over shorter intervals, comparable to the 
duration of Chandra observations of external galaxies. The re- 
sults of this analysis are presented in section |6] However, the 
effect of such time averaging may be more significant for tran- 
sients that, for long averaging times, will "accumulate" in the 
low-luminosity bins and will lead to a steepening of the XLF 
(IVoss & Gilfanovl2 007a). It is also an issue for M31, which was 
observed in several short observations distributed over the time 
span of a few years, and to a lesser extent for M81. This can 
potentially lead to significant distortions of the XLF, depend- 
ing on the average time and light curve properties of transients. 
However, as discussed in section |6^ it is not a significant factor 
in our particular case. 

It is also known that on e of the Galactic GCs, M l 5, contains 
two bright X-ray binaries dWhite & Angeli ni"200r). Transients 
hav e been detected in two other GCs, N GC 6440 (Heinke et al. 
12010 ) and Terzan 5 (iBordas et al.ll20ldl) . The net effect of source 
blending on the XLF of GCs was considered and shown to be 
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(erg s-1) 

Fig. 5. The incompleteness functions for individual galaxies 
from our sample (top panel for GC-LMXBs and middle panel 
for field LMXBs; the numbers 1, 2, 3, 4, 5, 6, and 7 by the 
curves refer to M31, Centaurus A, M81, Maffei 1, NGC 3379, 
NGC 4697, and NGC 4278, coiTespondingly). The bottom panel 
shows the combined curves for GC (solid line) and field (dotted 
line) sources. Lx is the unabsorbed X-ray luminosity in 0.5-8 
keV. 



negligible in IVoss & Gilfanovl (l2007ah . It is also discussed in 
section |6] 

4. Combined XLFs of GC and field LIVIXBs 

In combining the data from individual galaxies two effects need 
to be considered: correction for the incompleteness and removal 
of the contamination by background AGN. The former has to 
compensate for the fact that different sensitivities have been 
achieved for different galaxies, as well as for variations in the 
source detection sensitivity across the FOV in individual ob- 
servations. On the other hand, estimation of the contribution of 
CXB sources in the XLF has to take into account the effects of 
incompleteness that affect the detection of CXB sources as well. 

There are a few different weighting schemes for combin- 
ing the XLFs; here, we have used the one that produces the 
best signal-to-noise ratio under the assumption that XLFs of dif- 
ferent galaxies have the same shape (Eqs. ^ and (|6]l below). 
Obviously this assumption can only be verified to the accuracy 
allowed by the statistical quality of the individual XLFs, which 
is by a factor of a few less than the accuracy of the combined 



7 



Zhongli Zhang et al.: LMXB LPs in different environments 



XLF. Indirectly, this assumption is supported by not seeing large 
variations in the specific frequency of GC and field sources be- 
tween galaxies (Table|4] but see also section I6l2b . 

4.1. X-ray incompleteness correction 

We calculated the X-ray inc ompleteness function Kx(L) follow- 
ing the method of lVoss & Gilfanov. (.2006). For each pixel in the 
study area the point source detection sensitivity was calculated 
by inverting the detection method and using the actual local PSF 
and background (Voss & Gilfanov 2007a). Kx{L) is computed 
as the fraction of pixels, weighted by the assumed spatial dis- 
tribution of sources, in which a source with the given or higher 
luminosity would be detected. This computed quantity is the de- 
tection efficiency as a function of source luminosity. Because 
it depends on the spatial distribution of sources, we compute it 
separately for different source populations. For CXB sources we 
assume a flat spatial distribution. Field LMXBs are assumed to 
follo w the /T-band light , as reported in the 2MASS Large Galaxy 
Adas ('J arrett et al.l2003 ). For GC-LMXBs no weighting was ap- 
plied, but only pixels containing GCs were used in the calcula- 
tion. The incompleteness function for GC and field sources in 
individual galaxies are shown in the two upper panels in FigjS] 

The bottom panel in Fig|5] shows the combined incomplete- 
ness functions for GC and field sources computed by summing 
individual incompleteness functions weighted by the number of 
GC and stellar mass inside the study area of each galaxy. Since 
the distributions of GCs and stellar mass do not differ strongly, 
one may expect these two functions to be nearly identical. This 
is in fact the case throughout most of the luminosity range. The 
two curves diverge near ~ 5 x 10^"' erg s"' because of the dif- 
fere nt areas used to s tudy GC and field sources in Centaurus A 
(see IVoss et al.ll2009l for details). The difference below ~ 10^^ 
erg s^' is caused by only using GC sources for the MiUcy Way. 

4.2. LF of GC-LMXBs 

There are 185 GC-LMXBs in total in our sample. To avoid un- 
certainties from the highly incomplete low-luminosity end, we 
adopted a completeness threshold of 0.5 and used the curves 
shown in the upper panel of Fig |5] for each galaxy to determine 
the corresponding luminosity limit. Sources below these limits 
were excluded from the XLF construction. This procedure ex- 
cluded source numbers 13, 39, 98, 107, 108, and 109 from Table 
|6] The XLF value in the j-th luminosity bin centered at Lj and 
having a width of A log(L j) was computed according to the equa- 
tion: 



dN 



[dlogiLj) 



GC 



1 



V 



Alo 



eAL, 



1 



(if) 



(2) 



where N'~'£ is the effective number of GCs involved in the cal- 
culation of the XLF at the given luminosity Lx, corrected for 
optical and X-ray incompleteness. It depends on X-ray luminos- 
ity because it accounts for the X-ray incompleteness. 



K 

GC op' 



GC,X 



(3) 



where A^^*- is the number of observed GCs in the study field of 
the k-th galaxy, K^^{L) is the X-ray incompleteness function for 
GC sources in the k-th galaxy, the K^p,'f and K^p,i^ are optical 



completeness factors, described in section l331 and listed in Table 
[3] The thus computed XLF is normalized per GC. 

The effective X-ray incompleteness of the GC data can be 
defined as 



(4) 



This is the quantity plotted in Fig |5] (individual K{L) not clipped 
at the incompleteness level of 0.5 when plotting the figure). 

The factor K^^/^ /K^^, in eq.Q accounts for the incomplete- 
ness of the optical GC data. The denominator in this expression, 
Kgpi''^, is rather poorly constrained by our data (Table |3) and is 

consistent with K^p,'^ - K^p, within the measurement uncertain- 
ties. In fact, given the amplitude of the uncertainties, using the 
best-fit values would introduce additional noise into the obtained 
XLF. We therefore assumed that K^p,'^ - K^^, in the further cal- 
culations. With this assumption the incompleteness of the optical 
data cancels out (see the discussion in the beginning of the last 
paragraph in section l33t . The impact of this assumption on the 
final LF of GC-LMXBs is investigated in section 

In the case of GC sources, the contamination by background 
AGN is insignificant so was ignored. Indeed, the predicted to- 
tal number of random matches between resolved CXB sources 
and GC positions with the given search radii from Table |3] is 

1.1. The final XLF of GC sources is shown in Fig|6] The 
incompleteness-corrected number of GC sources with luminos- 
ity exceeding 10^''' erg s"' is ^ 244. 



4.3. LF of fieid LMXBs 

We have only considered field sources in elliptical galaxies and 
the bulges of spiral galaxies, in order to minimize the contam- 
ination by HMXBs. M81 was not included due to the rela- 
tively small size of its bulge. Thus, the bulge regions of M31, 
Maffei 1, NGC 3379, NGC 4278, and NGC 4697 were com- 
bined with the Centaurus A observations. For the last, we ex- 
cluded the jet and radio lobe r egions that have s mall-scale struc- 
tures in the diffuse emission (I Voss et al.ll2009h . In M31 we ex- 
cluded the sources located in the central 1 ' of the galaxy since 
they have been demonst rated as very likely dynamically formed 
(IVoss & Gilfanov! I2OO7 a) . A separate XLF was constructed for 
these 36 sources, as discussed in the next section. In Maffei 1, 
NGC 3379, NGC 4278, and NGC 4697 the central 10" were ex- 
cluded. These regions are affected by source confusion, which 
makes accurate luminosity estimates difficult. We then followed 
a procedure similar to the GC sources by applying a luminos- 
ity threshold corresponding to K(L) = 0.5 in each galaxy. With 
these selection criteria we obtained 496 sources above lO-'^ erg 
s"', of which ~177 are predicted to be CXB sources. 

In order to correctly subtract CXB contribution one has to 
take into account the difference in the incompleteness functions 
for CXB sources and LMXBs: 



dN 



dlogiLj) 



Alog( 



1 ( 1 

y I 



j 

JLe&Li 



A7tD\ 



dNcxB K^'^'il^) ^ 



dL K, 



LMXB 



(L) 



-dL 



(5) 



where Dk is the distance to the k-\h galaxy, K," {L) the com- 
bined incompleteness function for LMXBs computed similar 
to eq.Q, and dN'^^^/dS the log(A^) - log(5) distribution for 
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Fig. 6. The combined XLFs of LMXBs in GC (filled circles) and 
in the field (open circles). is the unabsorbed X-ray luminosity 
in 0.5-8 keV. The contribution of CXB sources was subtracted 
and the incompleteness correction was applied. The field XLF is 
normalized to stellar mass of 10'" Mq. The GC XLF is normal- 
ized to have the same number of sources as the field XLF above 
10^*^ ergs"'. 

Table 5. Comparison of XLFs of different populations. 



lO'- 



10^' 





89+""* 




Q 7, +0.13 
"■'-'-0.12 


Rf 






1 6T°-^^ 

l.u/_Q27 


Rc 






1 Q 1+0.97 
^•'^-0.58 


P(Rf < Rgc) 


<10-' (>5o-) 


3.6 


• 10-'' (3.6or) 


P{Rc < Rgc) 


8.3 • 10-3 (2.60-) 


4.6 


• 10-3 (2.80-) 



Note 



- Rgc-, Rf, and Rc are ratios of the number of faint sources 
< L, < 10^^ erg s-') to the number of bright sources (L, > 10^' 



the probabihty that the luminosity distributions of corresponding popu- 
lations are drawn from the same mean (see text for details). 



the CXB sources. In practice we implemented this by adding 
a large number (^lO"') of fake sources with small negative 
weights to each galaxy's source list. This accounted for the CXB 
log(A^) - log(5' ) distribution and incompleteness function of the 
galaxy. The sum of these weights for each galaxy equals the pre- 
dicted number of CXB sources in this galaxy. These "enhanced" 
source lists were used to produce the combined XLF according 
to eq.©: 



dN 



[dlogiLj) 



LMXB 



1 1 

v-LMXi 



LMXB^jJc^ 
\ 



■ z 



cxb 



A, 



LMXB rj cxb 



(6) 



The final XLF for field sources normalized to unit stellar mass 
is shown in Fig|6l The total stellar mass involved in this calcula- 
tion is 1.82x10"Mo. The specific frequency of LMXBs above 
10^^(10^^) erg S-' is 25.7(9. 6) per 10' ' ^Mo, w hich is consistent 
with the average values from lGilfanovl (12004 - 33.9(14.3). 



5. Results 

The background-subtracted and incompleteness-corrected XLFs 
of the GC and field LMXBs are shown in Figs. |6] and [T] The 
XLF of field sources is normalized to a stellar mass of 10'° Mq. 
The XLF of GC sources is noimalized to the same number of 
sources above lO-''^ erg/s as the field XLF. It is obvious from the 
plot that the two luminosity distributions have different shapes. 
Although they differ across the entire luminosity range, the most 
evident difference is at lower luminosities, below log(L) ~ 37. 
Both XLFs change their slope between log(L) ~ 37 - 38. Due to 
their rather complicated shapes we did not attempt to do global 
fits with analytical functions. Instead, we perform power-law 
fits to the high- and low-luminosity ends. We did maximum- 
likelihood fits to the background-subtracted XLFs. To account 
for the incompleteness, we multiplied the model by the respec- 
tive incompleteness function. In the log(L) > 38 range, we ob- 
tained diff^erential slopes of 1.70+°-^|| and 2.06+°-^2 ^j. qq and 
field sources coiTespondingly. At the faint end, log(L) < 37 the 
slopes are: 0.68^02' and l.n;^^}^ respectively. The slopes of the 
field sources are broadl y consistent with the parameters of aver- 
age LMXB XLFs froml Gilfanovl (12004 . 

Differences in the incompleteness curves and in the CXB 
contribution render direct application of the K-S test to com- 
pare these two XLFs impossible. We have therefore considered 
the ratio of the number of faint to bright sources in order to 
assess the statistical significance of the difference between the 
two XLFs, the same method as used in Voss et al. (2009). For 
each population we computed the ratio R - A^faint/Minght, with 
the boundary between faint and bright set to 10^^ erg s-'. We 
ran Monte-Caiio simulations to calculate statistical eiTors and 
the significance of our resul t s. The details of these calculations 
are described in IVoss et al.l (l2009l) . For each XLF we did 10^ 
Monte-Carlo runs. The results are listed in Table |5] The lines 
marked "P" give the probability of obtaining the observed values 
of R due to statistical fluctuations, while their mean (true) val- 
ues obey the relation given in parenthesis. These numbers can be 
interpreted as the probability that the corresponding luminosity 
distributions are drawn from the same mean. The upper limit of 
< 10-^ in the left column means that no such realizations were 
detected in 10^ Monte-Carlo runs. 

These calculations show that the GC and field XLFs differ 
at a confidence level of < 10"^, which corresponds to a signif- 
icance of > 5cr. To investigate the robustness of this conclu- 
sion we have also used a more restricted luminosity range of 
\og{Lx) > 36, where the incompleteness functions vary less and 
different galaxies from our sample make more uniform contri- 
butions. In this case the confidence level decreases to 3.4 ■ lO-"* 
(3.6cr), but the conclusion still holds. 



6. Discussion 

Our analysis has revealed a significant difference in the luminos- 
ity distributions of field and GC sources. A s imilar result was 
previously reported for a few nearby gal axies (IVoss & Gilfanovl 
I2007al: IVoss et al.ll20()9l: iKim et al .''200^. some of which are in- 
cluded in our sample. However, limited numbers of sources in 
these galaxies (even in Centaurus A) did not permit any de- 
tailed investigation of the luminosity distributions, but merely 
suggested a deficit of faint sources among dynamically formed 
LMXBs. Given the sensitivity expected in a typical Chandra 
observation and distances to the closest massive galaxies, the 
main limitation of these studies - insufficient number of sources 
cannot be lifted when investigating a single galaxy. This moti- 



9 



Zhongli Zhang et al.: LMXB LPs in different environments 



10.0 



-5 1.0 




Fig. 7. The combined XLFs of LMXBs in different environments 
plotted in the cumulative form. Lx is the unabsorbed X-ray lumi- 
nosity in 0.5-8 keV. The contribution of CXB sources was sub- 
tracted and the incompleteness correction was applied. The field 
XLF (solid) is normalized to the stellar mass of 10'" Mq. The 
normalizations of GC (dashed) and M31 nucleus (dash-dotted) 
XLFs are arbitrary. The shaded areas around the curves show 1 cr 
statistical uncertainty. 



vated us to combine data for several nearby galaxies. Thus, we 
have produced average XLFs of different populations of LMXBs 
with far better statistical quality than achieved before. With these 
XLFs we confirm the general conclusion of previous studies that 
the fraction of faint sources in GCs is a factor of ~ 4 smaller 
than sources in the field. Moreover, the overall behavior is more 
complex and cannot be described in terms of a simple roll-over 
of the XLF at low luminosities. Instead, luminosity distributions 
of GC and field sources differ throughout the entire luminosity 
range. 

It has long been debated whether the entire population of 
LMXBs in galaxies, including those in the field, was formed dy- 
namically i n GCs (White et al. 2002; Kundu et al. 2002, 2007; 
Ilrwinll2005t Duett 2005; Humphrey & Buote 2008|)- In this sce- 
nario it is suggested that either field LMXBs are systems ex- 
pelled from their GCs, or the host cluster itself was destroyed 
leaving the LMXB to join the field population. The significantly 
different luminosity distributions of LMXBs residing in GCs and 
in the field argues against this scenario. As is well known, the 
mass transfer rate in a Roche lobe-filling system is defined by 
its orbital period, mass ratio, and the evolutionary stage of the 
donor star. The distributions of these parameters are obviously 
different in the population of primordially formed binaries and in 
the dynamically formed systems in GCs. It is therefore not sur- 
prising that these populations have different mass transfer rates 
distributions and, correspondingly, different XLFsQ 



^ As a caveat, we mention that the cross-section for the dynamical in- 
teraction and, therefore, the probabihty for a binary to be ejected from 
a GC, is a strong function of the orbital period of the binary. Therefore, 
in the GC-scenario one may also expect different orbital period distri- 
bution of the binaries retained in GCs and those ejected into the field. 
This may lead to different XLFs of the two sub populations, if ejec- 
tion of binaries is the main mechanism of populating the field LMXBs, 
rather than GC destruction. Detailed populations synthesis calculations 
are required in order to see if this may explain the observed XLFs. 
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Fig. 8. The XLFs of sources in the inner 1 ' of M3 1 (open circles) 
and of GC sources (filled circles). L^ is the unabsorbed X-ray 
luminosity in 0.5-8 keV. The normalization of the GC XLF is 
arbitrary. No sources are detected in the three highest luminosity 
bins of the M31 XLF 



6. 7. Sources in the nucleus of M31 



IVoss & Gilfanovl (l2007ah demonstrate that most sources in the 
inner ~ 1 ' of M3 1 are very likely to have formed dynamically, 
similar to the sources in GCs. In particular, they find that their 
spatial distribution follows the law, in contrast to the X-ray 
sources outside this region, where the density is proportional 
to the stellar density. We updated the LMXB list in this region 
using an exposure appro ximately twice of the one presented in 
IVoss & Gilfanovl (l2007a). We have detected three new sources, 
bringing the total number to 36 (The increase in the number of 
sources for a log(A^) - log(5' ) distribution with the slope of - 1 
would be ~13). We excluded one source that coincided with a 
GC and computed the luminosity distribution of the detected 
sources, performing incompleteness correction and CXB sub- 
traction as described before. The resulting XLF is shown in Fig. 
|8]along with the XLF of GC sources. It is obvious from the plot 
that the two distributions have similar shapes at log(Lx) ~ 37 
but differ at the bright end, with the XLF of the sources in the 
M31 nucleus having a deficit of bright sources. To test the sta- 
tistical significance of this conclusion we ran the same tests as 
we did to compare XLFs of GC and field LMXBs in Sec|5] We 
found that the LF of LMXBs in the M31 nucleus differs from 
the GC XLF with a significance ai2.6-2.8cr (Table |5]i. In other 
words, one should expect 26 sources with log(Lx) > 37 in the 
nucleus of M31, assuming that both distributions have the same 
shape and using the number of faint sources for normalization. 
The observed number of bright sources is 11, which is x 3cr less 
than expected. Both calculations give similar results, confirming 
the marginal significance of our conclusion. 

Although both populations were formed dynamically, there 
is an important difference between stellar environments in GCs 
and galactic nuclei: stellar velocities in the latter are about 
~ 10-20 times higher. This l eads to diffe rent formation chan- 
nels in GCs and galactic nuclei (IVoss & Gilfa nov 2007b, and ref- 
erences therein). Calculations of IVoss & Gilfanov. (■2007b ) sug- 
gest that in the high-velocity environment of the M3 1 nucleus 
the main formation channel for X-ray binaries may be tidal cap- 
tures of compact objects by low-mass stars, producing short or- 
bital period binaries. In GCs, in contrast, LMXBs are predomi- 
nantly formed in exchange reactions and collisions of neutron 
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Fig. 9. The maximum possible effect of the optical incomplete- 
ness on the GC LMXB XLF. The solid symbols show the XLF 
computed assuming 

/r^CA- ^ J.CC (qjjj. default version), open 

circles - assuming K^^/^ — 1 (the maximum possible correc- 
tion). To emphasize the effect on the shape, rather than overall 
normalization, the XLF are normalized to the total number of 
detected GCs. See sections 1331 14.21 and l6.2l for details. 



stars with evolved stars. Obviously, this difference will affect 
the distributions of binary systems over the mass accretion rate. 
Detailed population synthesis calculations are required to un- 
derstand what this effect may be and to interpret the observed 
luminosity distributions in a more quantitative and meaningful 
manner. 

6.2. Caveats 

Several caveats regarding possible systematic effects are in or- 
der. The distance uncertainties, in the limit of several galaxies, 
will smooth out the luminosity distributions. For a smaller sam- 
ple (which is the case for the current study), there could be a 
non trivial effect on the computed XLFs. However, the distances 
to the galaxies in our sample are fairly well known, with an ac- 
curacy of ~ 5 - 15% (Table [TJ. This translates to ~ 10 - 30% 
uncertainty in the luminosity and ~ 4 - 11 % in its logarithm. 
This is a factor of ~ 4 - 12 smaller that the bin width used in the 
XLF calculations. Thus it should not affect the measured XLFs 
in any significant way. 

We used a rather inhomogeneous set of the GC lists, having 
varying degree of completeness. The procedure of correction for 
incompleteness of the optical data is described in section lTD and 
14.21 Its accurate implementation, however, was hindered by the 
large statistical uncertainties of the completeness fraction of the 
GCs hosting an X-ray source, K^^,'^ . As results of section [33] 

were consistent with K^p,'^ - K^^, (the latter much better con- 
strained), we assumed that this relation holds for all galaxies. 
This could be the case, for example, if the probability of finding 
an LMXB did not depend on the optical luminos ity of the GC. 
Howe ver, this is known not to be the case (e.g., ISivakoff et alJ 
l2007h (see Fig.|4]i. As is obvious from Eq.©, the optical incom- 
pleteness would have the strongest effect on the XLF if for all 
galaxies K^^,'^ = 1 (i.e. if all LMXBs were located in the bright- 
est GCs and not subject to optical incompleteness at all, which is 
not true either). To illustrate its amplitude, we show in Fig |9] the 
XLFs computed in these two limiting cases. As is obvious from 
the plot, the XLF does not change by more than ~ 20 - 30% in 



the two lowest luminosity bins. We emphasize that the example 
shown in the plot illustrates the maximum possible effect of the 
optical incompleteness, the real effect being smaller. 

Combining XLFs necessarily involves an assumption regard- 
ing the similarity of their shapes in individual galaxies. Although 
we did not detect statistically significant differences between dif- 
ferent galaxies, this assumption cannot be verified directly at the 
same level of accuracy as provided by the output average XLF. 
On the other hand, we do detect marginally significant varia- 
tions in the specific frequency of X-ray sources in GCs between 
galaxies, although these may be related, at least in part, to the in- 
completeness of the GC lists in more distant galaxies. However, 
if they are real, they may be accompanied by variations in the 
XLF shapes. The effect of such variations may be further am- 
plified by the fact that data for different luminosities come from 
different galaxies. The low-luminosity domain, log(L) < 36, is 
covered exclusively by the nearby M3 1 and Milky Way, whereas 
the bright end is dominated by sources located in more mas- 
sive but more distant galaxies, such as Centaurus A and NGC 
4697. This is another unavoidable limitation, as bright sources, 
although more easy to detect, are less frequent, and it takes a 
bigger galaxy to have them in large numbers. On the other hand, 
bigger galaxies are more distant and the sensitivity achieved in 
a typical Chandra observation is lower. Conversely the nearby 
galaxies, where fainter sources can be studied, tend to be less 
massive and contain fewer blight sources. Luckily, the 800 ksec 
Very Large Chandra program on Centaurus A and relatively 
good coverage of M3 1 allowed us to bridge faint, intermediate, 
and bright luminosity ranges. 

About half of the GCs with X-ray sensitivity in the lowest 
luminosity domain, log(L) ^ 36 are located in the Milky Way. 
The flux determination of the latter, based on the averaging of 
the ASM light curves, may be subject to systematic effects. It 
is different from those affecting Chandra galaxies data as dis- 
cussed in section IT?! Piimaiily, this is due to the long integration 
times of ASM light curves. To investigate its effect on the XLF, 
we divided the ASM light curves into 100 sub intervals with a 
duration of 50 days each (comparable to the integration time of 
the longest Chandra observations) and recalculated the GC XLF 
100 times, each time using the data from different sub intervals 
to compute ASM fluxes for the Milky Way sources. The range 
of obtained XLF values is shown by shaded area in Fig.[TO] As 
is obvious from the plot, the long integration time of ASM data 
does not affect the GC XLF significantly. 

If the time span of observations is longer than the typi- 
cal time scale of transient sources, averaging of their luminos- 
ity can also modi fy the shape of the lu minosity distribution, 
making it steeper ( Voss & Gilfanov 2007a). The typical decay 
time scales of the Galactic transient sources are in the ~weeks- 
months range. Thus, for the Milky Way GCs, this issue is ad- 
dressed by the above exercise with the ASM light curves, and 
Fig [To] demonstrates that averaging of transients does not result 
in significant modifications of XLF, given its statistical quality. 
This issue is also relevant for the multiple Chandra observations 
of M3 1 and, to a less extent, M8 1 . Indeed, the Chandra image 
of the bulge of M3 1 was obtained by combining more than 40 
short (~5 ks) observations. As transients are bright only in a 
few observations and dim in many others, they will tend to ac- 
cu mulate in the low-lumino sity bins, making the XLF steeper. 
In I Voss & Gilfanov! (l2007ah 28 transients were reported, two of 
which (Src. 22 and Src. 35 in Table |6) are in our GC-LMXB 
source list and 21 are in the field source list. We recomputed the 
luminosity distributions excluding these sources and did not find 
any significant changes (Fig. [TTI) . The results of the statistical 
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Fig. 10. The combined XLF of GC-LMXBs. Lx is the unab- 
sorbed X-ray luminosity in 0.5-8 keV. The XLF uncertainty 
shown by the shaded regions is due to the variability of GC- 
LMXBs in the Milky Way. 
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Fig. 11. The combined XLFs of field sources with (open circles) 
and without (filled circles) transient sources in M31. Lx is the 
unabsorbed X-ray luminosity in 0.5-8 keV. Due to multiple short 
observations averaging the luminosities of transients (which are 
bright in only a few observations and faint in many others) the 
number of faint sources is artificially increased, making the XLF 
appear steeper See text for details. 



tests reported in Table |5]are not changed significantly either: the 
probability P(Rf < Rcc) for the full luminosity range remains 
< 10"^ (although this may be affected somewhat by the incom- 
pleteness of the transient list at the faint end) for sources with 
log(Lx) > 36 it changes from 3.6 ■ lO"'* to 2 ■ 10"^ (from 3.6cr to 
3.1cr). This proves that the contamination by transient sources in 
M31 does not significantly bias our results. It is much less sig- 
nificant for other galaxies, as they were observed by Chandra in 
much fewer longer observations. Thus, averaging of persistent 
and transient sources does not lead to significant (as compared 
to statistical errors) distortions of the XLFs derived in this paper 
and does not affect our results in any significant way. This con- 
clusion should not be taken out of the context though. In a more 
general case, the effects discussed above may be important and 
may need a more elaborate treatment. 

Another factor that can modify the apparent XLF of GC 
LMXBs is the multiplicity of X-ray sources in GCs, which can 
affect both Milky Way data and Chandra data for external galax- 



ies. We use the Milky Way GCs to estimate its impact on XLF. 
One of the Milky Way GCs, M15, is known to contain two 
persistent LMXBs (X2127 -1-119-1 and X2127+119-2) with in- 
stantaneous luminosities of 9.5x10^^ erg s~' and 1.5x10-^^ er g 
s ' (converted to 0.5-8 keV band) (iWhite & Angelinill2001h . 
Obviously these two sources could not be resolved by ASM, 
which measured the long-term average luminosity of 4.05x10"'^ 
erg s"'. Similarly, they would not be resolved by Chandra in 
any of our external galaxies with the exception of M31, where 
it may be marginally possible. Two b right transient s have been 
de tected recently, in NGC 6440 (Hein ke eta"l]|2010h and Terzan 
5 (iBordas et al.ll20i 0) with the luminosities in the ~ lO^*" - 10^** 
erg/s range. Their effect on the "snapshot" XLF would depend 
on their unknown duty cycle. Assuming a duty cycle of ~ 0.5, 
which seems to be a very generous upper limit, the multiplicity 
fraction for the Mil ky Way GCs is -1/12- 2/12 ^ 8 - 16%. 
As demonstrated in lVoss & Gilfanovl (l2007a). the multiplicity at 
the level of ~ 10% does not modify the luminosity distribution 
in any significant way. We also checked to see how the GC XLF 
is affected if the ASM M15 source is replaced by two source 
with the luminosities determined by Chandra and the transient 
source in NGC 6440 at its brightest state is added to the sam- 
ple. The overall effect on the XLF is insignificant and the Rqc 
changed from 0.89+[5 [^ to 0.92+° }^, which is also neghgible. We 
conclude that unless the multiplicity is much higher in external 
galaxies, it does not affect our conclusions in any significant way. 



7. Summary 

The aim of this study was to produce accurate luminosity dis- 
tributions of LMXBs in different environments - dynamically 
formed systems in GCs, in the nucleus of M3 1 and field sources 
of presumably primordial origin - in order to facilitate their 
quantitative comparison and to provide input for verifying pop- 
ulation synthesis models. This goal required a broad luminos- 
ity coverage with a point source detection sensitivity reaching 
10^^ erg/s and, on the other hand, good sampling of the high- 
luminosity end, where the specific frequency of sources (per GC 
or per unit stellar mass) is low. As this combination of proper- 
ties cannot be achieved with a single galaxy, we combined the 
data from a number of galaxies. To this end, we assembled a 
sample of galaxies from the public Chandra archive which is 
best suited to our study. It included seven nearby galaxies (M31, 
Maffei 1, Centaurus A, M81, NGC 3379, NGC 4697, and NGC 
4278) and the Milky Way. We detected 185 X-ray sources in 
1615 GCs, 36 sources in the nucleus of M31, and 998 sources in 
the fields of galaxies (of which » 365 are expected to be back- 
ground AGN). These sources were used to produce the average 
luminosity distributions of different populations. In doing so we 
took special care to accurately subtract resolved CXB sources 
and correct for incompleteness effects. As a result, we produced 
XLFs of LMXBs with a statistical accuracy that far exceeds what 
has been achieved in previous studies. 

We demonstrate that, although the luminosity distributions 
of LMXBs in different environments are similar in a broad 
sense (e.g., when compared with XLF of HMXBs), their de- 
tailed shapes are different. Although the fraction of faint LMXBs 
(log(Lx) < 37) in GCs is ~ 4 times smaller than in the field, in 
agreement with a suggested effect found in previous studies, the 
difference in their XLFs cannot be described merely in terms of 
a roll over of the XLF of GC sources. Rather, the luminosity dis- 
tributions of these two populations of LMXBs appear to be dif- 
ferent throughout the entire luminosity range. This may present 
a challenge for the models suggesting that the entire LMXB pop- 
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ulation was formed dynamically in GCs and then expelled to the 
field due to kicks, dynamical interactions, or GC destruction. 

We also compare luminosity distributions of LMXBs in the 
nucleus of M31 (its inner 1') and in GCs. We find that although 
their shapes at the low-luminosity end are similar (and dififer- 
ent from the field sources), the M3 1 nuclear population appears 
to have far fewer luminous sources than GCs (and field popula- 
tion). For example, the most luminous source in the nucleus of 
M31 has the luminosity of 4.7 x 10^^ erg/s. If the XLFs were 
drawn from the same parent distribution, we would expect to 
see 1 1 sources above this luminosity, whereas we found none. 
Different estimates of the statistical significance of the differ- 
ence between the two XLFs give results in the ~ 2.5 - 3cr range. 
The difference between the XLFs is likely caused by the factor 
of ~ 10-20 difference in stellar velocities in GCs and galactic 
nuclei, which leads to different dynamical formation channels. 
However detailed population synthesis calculations are needed 
in order to understand the particular mechanisms responsible for 
forming the observed luminosity distributions. 
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Table 6. All the 185 LMXBs in GCs in our sample. 



Number 


Galaxy 


RA(J20()0) 


DEC(J2000) 


Luminosity 


Number 


Galaxy 


RA(J2000) 


DEC(J2000) 


Luminosity 


(1) 


(2) 


(3) 


(4) 


(5) 


(1) 


(2) 


(3) 


(4) 


(5) 


1 


MW 


+ 17:35:47.64 


-30:28:55.70 


0.14 


74 


Cen A 


+ 13:26:00.72 


-43:09:40.32 


46.48 


2 


MW 


+18:53:04.89 


-08:42:19.70 


1.14 


75 


Cen A 


+ 13:25:46.56 


-42:57:02.88 


55.91 


3 


MW 


+ 17:48:53.54 


-20:22:02.00 


1.31 


76 


Cen A 


+ 13:25:35.52 


-42:59:34.80 


57.54 


4 


MW 


+ 18:35:44.00 


-32:58:55.40 


1.68 


77 


Cen A 


+ 13:25:09.12 


-42:58:58.80 


58.10 


5 


MW 


+ 17:50:45.54 


-31:17:32.50 


1.76 


78 


Cen A 


+ 13:25:10.32 


-42:53:32.64 


69.64 


6 


MW 


+ 17:48:55.73 


-24:53:40.10 


1.78 


79 


Cen A 


+ 13:25:12.96 


-43:01:14.16 


79.97 


7 


MW 


+05:14:06.59 


-40:02:37.00 


2.34 


80 


Cen A 


+ 13:25:31.68 


-43:00:02.88 


87.42 


8 


MW 


+21:29:58.33 


+ 12:10:02.80 


4.05 


81 


Cen A 


+ 13:25:35.28 


-42:53:00.96 


97.94 


9 


MW 


+ 17:33:24.06 


-33:23:16.20 


5.00 


82 


Cen A 


+ 13:25:54.48 


-42:59:25.44 


105.59 


10 


MW 


+ 17:27:33.25 


-30:48:07.40 


6.79 


83 


Cen A 


+ 13:25:07.68 


-43:01:14.88 


199.09 


11 


MW 


+ 17:50:12.66 


-37:03:08.20 


14.70 


84 


Cen A 


+ 13:25:02.64 


-43:02:43.08 


255.53 


12 


MW 


+ 18:23:40.57 


-30:21:40.60 


66.12 


85 


Cen A 


+ 13:25:42.00 


-43:10:41.52 


315.01 


13 


M31 


+00:42:29.64 


+41:17:57.27 


0.04 


86 


M81 


+09:56:05.30 


+69:06:43.53 


2.01 


14 


M31 


+00:42:50.86 


+41:10:33.72 


0.41 


87 


M81 


+09:55:37.26 


+69:02:07.57 


2.36 


15 


M31 


+00:43:14.65 


+41:25:13.32 


0.84 


88 


M81 


+09:55:51.97 


+69:07:39.18 


4.83 


16 


M31 


+00:42:27.43 


+40:59:35.63 


1.05 


89 


M81 


+09:55:22.05 


+69:05:18.93 


6.68 


17 


M31 


+00:42:34.40 


+40:57:09.31 


1.07 


90 


M81 


+09:55:54.93 


+69:00:56.03 


35.33 


18 


M31 


+00:43:15.48 


+41:11:25.69 


1.17 


91 


M81 


+09:55:47.00 


+69:05:51.09 


67.26 


19 


M31 


+00:42:40.60 


+41:10:33.60 


1.39 


92 


M81 


+09:55:58.54 


+69:05:26.04 


70.40 


20 


M31 


+00:42:25.04 


+40:57:18.78 


2.09 


93 


M81 


+09:55:49.80 


+69:05:31.93 


434.56 


21 


M31 


+00:42:41.43 


+41:15:23.71 


2.53 


94 


Maffei 1 


+02:36:37.26 


+59:39:15.50 


5.41 


22 


M31 


+00:42:47.81 


+41:11:13.66 


2.56 


95 


Maffei 1 


+02:36:30.84 


+59:39:34.70 


9.11 


23 


M31 


+00:43:07.51 


+41:20:19.44 


3.24 


96 


Maffei 1 


+02:36:26.03 


+59:39:06.91 


16.22 


24 


M31 


+00:42:33.10 


+41:03:29.86 


4.23 


97 


Maffei 1 


+02:36:36.50 


+59:38:42.03 


42.07 


25 


M31 


+00:42:09.5 1 


+41:17:45.42 


9.31 


98 


N3379 


+ 10:47:50.47 


+ 12:34:23.11 


4.14 


26 


M31 


+00:42:26.05 


+41:19:14.81 


9.94 


99 


N3379 


+ 10:47:51.57 


+ 12:35:36.01 


14.18 


27 


M31 


+00:42:12.17 


+41:17:58.62 


11.51 


100 


N3379 


+ 10:47:54.20 


+ 12:35:29.49 


39.83 


28 


M31 


+00:43:03.31 


+41:21:21.60 


12.02 


101 


N3379 


+ 10:47:50.47 


+ 12:34:36.94 


53.20 


29 


M31 


+00:42:31.25 


+41:19:38.78 


18.50 


102 


N3379 


+ 10:47:50.33 


+ 12:35:06.59 


58.00 


30 


M31 


+00:43:02.93 


+41:15:22.47 


22.50 


103 


N3379 


+ 10:47:51.08 


+ 12:35:49.25 


87.35 


31 


M31 


+00:43:03.86 


+41:18:04.79 


28.37 


104 


N3379 


+ 10:47:52.77 


+ 12:35:08.58 


242.70 


32 


M31 


+00:42:59.86 


+41:16:05.64 


33.75 


105 


N3379 


+ 10:47:50.19 


+ 12:34:55.34 


333.36 


33 


M31 


+00:42:59.65 


+41:19:19.18 


34.19 


106 


N3379 


+ 10:47:52.65 


+12:33:38.01 


680.58 


34 


M31 


+00:42:18.64 


+41:14:01.74 


36.23 


107 


N4697 


+ 12:48:32.94 


-05:47:04.04 


7.23 


35 


M31 


+00:43:14.31 


+41:07:19.68 


46.14 


108 


N4697 


+ 12:48:33.63 


-05:48:49.20 


8.22 


36 


M31 


+00:43:37.29 


+41:14:43.63 


47.94 


109 


N4697 


+ 12:48:29.13 


-05:48:22.15 


8.26 


37 


M31 


+00:43:10.61 


+41:14:51.24 


77.93 


110 


N4697 


+ 12:48:34.64 


-05:47:27.55 


9.39 


38 


M31 


+00:42:15.84 


+41:01:14.32 


123.26 


111 


N4697 


+ 12:48:26.52 


-05:47:24.91 


9.53 


39 


Cen A 


+ 13:25:41.76 


-42:57:00.00 


0.90 


112 


N4697 


+ 12:48:35.80 


-05:47:41.90 


10.13 


40 


Cen A 


+13:25:1 1.04 


-43:01:31.80 


1.79 


113 


N4697 


+ 12:48:37.60 


-05:47:49.79 


10.26 


41 


Cen A 


+ 13:25:29.28 


-42:57:46.80 


1.90 


114 


N4697 


+ 12:48:34.68 


-05:48:14.82 


1 1 .08 


42 


Cen A 


+ 13:25:14.88 


-43:00:48.96 


1.92 


115 


N4697 


+ 12:48:37.16 


-05:48:30.34 


1 1.76 


43 


Cen A 


+ 13:25:42.00 


-43:03:19.44 


2.15 


116 


N4697 


+ 12:48:28.04 


-05:48:32.66 


13.69 


44 


Cen A 


+ 13:25:58.32 


-43:08:06.72 


3.15 


117 


N4697 


+ 12:48:40.86 


-05:48:23.12 


14.88 


45 


Cen A 


+ 13:25:43.20 


-42:58:37.20 


3.57 


118 


N4697 


+ 12:48:41.66 


-05:48:47.04 


15.00 


46 


Cen A 


+ 13:25:35.28 


-43:05:29.40 


4.45 


119 


N4697 


+ 12:48:37.71 


-05:47:29.32 


16.78 


47 


Cen A 


+ 13:24:49.20 


-43:05:12.12 


4.81 


120 


N4697 


+ 12:48:35.80 


-05:46:40.69 


17.44 


48 


Cen A 


+ 13:25:27.60 


-43:05:24.72 


5.19 


121 


N4697 


+ 12:48:31.84 


-05:48:38.70 


23.44 


49 


Cen A 


+ 13:25:32.40 


-43:04:40.44 


6.25 


122 


N4697 


+ 12:48:36.97 


-05:47:32.61 


23.97 


50 


Cen A 


+13:24:58.08 


-42:56:10.32 


6.84 


123 


N4697 


+ 12:48:33.95 


-05:48:34.46 


25.58 


51 


Cen A 


+13:25:14.16 


-43:02:42.72 


7.01 


124 


N4697 


+ 12:48:26.16 


-05:47:29.50 


26.29 


52 


Cen A 


+13:25:22.08 


-43:02:45.24 


7.68 


125 


N4697 


+ 12:48:36.95 


-05:48:10.80 


30.36 


53 


Cen A 


+ 13:25:30.24 


-42:59:34.80 


7.83 


126 


N4697 


+ 12:48:33.19 


-05:49:12.85 


40.98 


54 


Cen A 


+ 13:25:32.88 


-42:56:24.36 


8.54 


127 


N4697 


+ 12:48:37.87 


-05:46:52.81 


42.12 


55 


Cen A 


+ 13:24:50.40 


-43:04:50.88 


9.16 


128 


N4697 


+ 12:48:40.92 


-05:47:31.44 


42.31 


56 


Cen A 


+ 13:25:38.40 


-42:57:19.80 


10.12 


129 


N4697 


+ 12:48:31.05 


-05:48:28.66 


46.15 


57 


Cen A 


+ 13:25:12.00 


-42:57:12.96 


10.27 


130 


N4697 


+ 12:48:41.50 


-05:47:37.25 


46.82 


58 


Cen A 


+ 13:26:07.68 


-42:52:01.56 


10.40 


131 


N4697 


+ 12:48:36.10 


-05:48:33.61 


60.54 


59 


Cen A 


+ 13:26:05.28 


-42:56:32.64 


10.71 


132 


N4697 


+ 12:48:46.55 


-05:48:12.02 


75.50 


60 


Cen A 


+ 13:26:10.56 


-42:53:43.08 


11.14 


133 


N4697 


+ 12:48:38.67 


-05:47:46.88 


91.24 


61 


Cen A 


+ 13:25:05.76 


-43:10:30.36 


11.57 


134 


N4697 


+ 12:48:35.95 


-05:45:51.79 


91.79 


62 


Cen A 


+ 13:25:28.08 


-43:04:01.92 


13.71 


135 


N4697 


+ 12:48:31.73 


-05:48:46.73 


97.43 


63 


Cen A 


+ 13:25:03.12 


-42:56:24.72 


13.84 


136 


N4697 


+ 12:48:36.97 


-05:48:01.04 


110.29 


64 


Cen A 


+ 13:25:32.88 


-43:04:28.92 


15.85 


137 


N4697 


+ 12:48:32.65 


-05:48:51.11 


125.21 


65 


Cen A 


+ 13:25:52.80 


-43:05:46.32 


20.97 


138 


N4697 


+ 12:48:35.97 


-05:47:56.56 


150.11 


66 


Cen A 


+ 13:25:05.04 


-43:01:32.88 


22.37 


139 


N4697 


+ 12:48:39.35 


-05:47:30.48 


168.46 


67 


Cen A 


+ 13:25:39.84 


-43:05:01.68 


23.62 


140 


N4697 


+ 12:48:36.72 


-05:47:31.89 


178.73 


68 


Cen A 


+ 13:25:32.40 


-42:58:49.80 


23.69 


141 


N4697 


+ 12:48:37.51 


-05:47:43.40 


192.56 


69 


Cen A 


+ 13:25:18.48 


-43:01:15.96 


24.05 


142 


N4697 


+ 12:48:27.03 


-05:49:25.25 


206.76 


70 


Cen A 


+ 13:25:10.56 


-43:06:24.12 


27.47 


143 


N4697 


+ 12:48:30.83 


-05:48:36.93 


308.63 


71 


Cen A 


+ 13:26:19.68 


-43:03:19.08 


30.51 


144 


N4697 


+ 12:48:33.21 


-05:47:41.90 


455.00 


72 


Cen A 


+13:25:19.92 


-43:03:09.72 


31.62 


145 


N4697 


+ 12:48:39.32 


-05:48:07.22 


474.28 


73 


Cen A 


+13:25:10.32 


-42:55:09.48 


37.62 


146 


N4278 


+ 12:20:04.55 


+29:18:19.33 


7.97 
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Table 6. continued. 



Number 


Galaxy 


RA(J2000) 


DEC(J2000) 


Luminosity 


Number 


Galaxy 


RA(J2000) 


DEC(J2000) 


Luminosity 


(1) 


(2) 


(3) 




(4) 


(5) 


(1) 


(2) 


(3) 




(41 


(5) 


147 


N4278 


+ 12:20:00.39 


+29 


17:46.37 


8.53 


167 


N4278 


+ 12:20:07.16 


+29 


17:38.74 


50.30 


148 


N4278 


+ 12:20:04.70 


+29 


16:07.46 


9.19 


168 


N4278 


+ 12:20:04.53 


+29 


16:12.19 


50.37 


149 


N4278 


+ 12:20:02.98 


+29 


18:14.97 


9.78 


169 


N4278 


+ 12:20:00.32 


+29 


17:05.11 


52.01 


150 


N4278 


+ 12:20:00.37 


+29 


17:22.08 


11.84 


170 


N4278 


+ 12:20:03.77 


+29 


1 fi-09 66 


58.85 


151 


N4278 


+ 12:20:05.24 


+29 


16:01.51 


12.03 


171 


N4278 


+ 12:20:08.04 


+29 


16:42.13 


61.23 


152 


N4278 


+ 12:20:04.87 


+29 


16:01.73 


15.23 


172 


N4278 


+ 12:20:09.15 


+29 


17:57.95 


66.65 


153 


N4278 


+ 12:20:03.54 


+29 


16:17.50 


17.76 


173 


N4278 


+ 12:20:04.11 


+29 


16: 15.34 


66.84 


154 


N4278 


+ 12:20:02.49 


+29 


16:24.65 


18.28 


174 


N4278 


+ 12:20:08.07 


+29 


16:43.61 


71.32 


155 


N4278 


+ 12:20:05.89 


+29 


18:21.35 


18.72 


175 


N4278 


+ 12:20:08.85 


+29 


17:28.92 


90.19 


156 


N4278 


+12:20:01.08 


+29 


17:23.52 


21.70 


176 


N4278 


+ 12:20:08.39 


+29 


17:16.85 


113.49 


157 


N4278 


+ 12:20:04.59 


+29 


16:15.51 


22.01 


177 


N4278 


+ 12:20:05.70 


+29 


16:49.98 


119.54 


158 


N4278 


+ 12:20:09.95 


+29 


17:40.59 


24.65 


178 


N4278 


+ 12:20:08.15 


+29 


17:16.97 


125.36 


159 


N4278 


+ 12:20:05.07 


+29 


17:15.46 


26.87 


179 


N4278 


+ 12:20:05.95 


+29 


17:08.94 


138.53 


160 


N4278 


+ 12:20:06.33 


+29 


17:10.05 


28.26 


180 


N4278 


+ 12:20:07.71 


+29 


16:44.05 


144.98 


161 


N4278 


+ 12:20:08.14 


+29 


16:59.83 


28.32 


181 


N4278 


+ 12:20:06.82 


+29 


16:36.65 


145.01 


162 


N4278 


+ 12:20:00.28 


+29 


18:12.18 


29.04 


182 


N4278 


+ 12:20:05.24 


+29 


16:39.97 


265.89 


163 


N4278 


+ 12:20:01.85 


+29 


17:58.35 


30.38 


183 


N4278 


+ 12:20:04.23 


+29 


16:51.47 


269.32 


164 


N4278 


+ 12:20:05.24 


+29 


16:52.84 


37.40 


184 


N4278 


+ 12:20:03.44 


+29 


16:39.55 


292.63 


165 


N4278 


+12:20:02.00 


+29 


17:29.78 


46.63 


185 


N4278 


+ 12:20:07.76 


+29:17:20.46 


388.68 


166 


N4278 


+12:20:03.73 


+29 


16:29.81 


48.82 















(1) - the sequence number. (2) - the galaxy where source is detected. (3),(4) - Right ascension and declination of source. (5) - X-ray luminosity in 0.5-8 keV, 
in unit of 10^* erg s~'. 
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